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ABSTRACT 

We present the discovery and spectroscopic follow-up of a nearby late-type L dwarf (2M0614+3950), 
and two extremely wide very-low-mass binary systems (2M0525-7425AB and 2M1348-1344AB), re- 
sulting from our search for common proper motion pairs containing ultracool components in the Two 
Micron All Sky Survey (2MASS) and the Wide-field Infrared Survey Explorer (WISE) catalogs. The 
near-infrared spectrum of 2M0614+3950 indicates a spectral type L9 ± 1 object residing at a distance 
of 26.0 ±1.8 pc. The optical spectrum of 2M0525-7425A reveals an M3.0 ± 0.5 dwarf primary, accom- 
panied by a secondary previously classified as L2. The system has an angular separation of ~ 44", 
equivalent to ~ 2000 AU at distance of 46.0 ± 3.0 pc. Using optical and infrared spectra, respectively, 
we classify the components of 2M1348-1344AB as M4.5 ± 0.5 and T5.5 ± 1. The angular separation 
of ~ 68" is equivalent to ~ 1400 AU at a distance of 20.7 ± 1.4 pc. 2M1348-1344AB is one of only six 
very wide (separation > 1000 AU) systems containing late T dwarfs known to date. 
Subject headings: binaries:general - stars:low-mass, brown dwarf - stars: individual 

(2MASSJ06143818+3950357, 2M05254550- 7425263, 2M05253876- 7426008, 

2M13480721-1344321, 2M13480290-1344071 



1. INTRODUCTION 

The formation and evolution of objects at the low-mass 
end of the initial mass function is one of the most fun- 
damental and yet not fully understood issues in our cur- 
rent picture of how stars form and evolve. The existence 
of brown dwarfs (BDs; m < 0.08 M ) down to — 7Mj 
can be exp lained within the standard framework of star 
formation ( Low fc Lynden-Bell||1976[ ). By including ad- 
ditional ph ysics such as turbulence dPadoan fc NorcT 
lund| | 2004| , dynamical interactions (Boffin et al. |1998|" 
Bate||2009[). disk fragmentation (Stamatellos & Whit- 
worth"fp009j) , and/or photo-erosion from nearby bright 
stars flwhitworth fc Zinnecker 2004), one can create ob- 



jects with even lower masses (down to 3Mj). The rel- 
ative importance of these mechanisms in producing the 
Galactic substellar population is unclear. The key here is 
to confront predictions from the various scenarios with 
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observational results. One means of testing these the- 
ories is through the empirical characterization of very- 
low-mass ( VLM) multiple systems. As summarized in 
Burgasser (2007), these systems are observed to be less 
frequent, more tightly bound, and of higher mass ratios 
than their more massive counterparts. 

Wide binaries are of interest because their large sepa- 
rations (greater than 100 AU) and low binding energies 
provide direct constraints for formation models. For ex- 
ample, wide binaries are expected to be disrupted during 
the ejection process. This directly challenges the ejec- 
tion model for the formation of VLM stars and BDs 



( |Reipurth fc Clarke|200lj|Bate fc Bonnell|2005[ ). Several 
other models have been been successful in creating wide 



VLM bina ries, e.g. by fragmentation during the cluster 
formation (Bate 2009), in fragmenting disks (Stamatel- 
los & Whitworth 2009 ) , or by dynamical capture during 
star cluster dissolution p hase (Kouwenhoven et all 2010| 



and Kouwenhoven et al 



Mo eckel fc Clarke 2011 ) . The models by Bate 7 2009 



(20101 predict the frequency of 



binaries with separations > 1000AU to be a few percent 
compared with a fraction of ~ 10% found observation- 
ally for nearby (Lepine & B ongiorno 2007; Longhitano 



no 
9) 



& Binggcli 2010) and young ( jKraus fc Hillenbrand|200' 
stars. In star forming regions, the fraction of wide bi 
naries appears to decrease with decreasing stellar mass 
( Kraus & Hillenbrand 2009 ) . While the number of known 
wide binaries containing ultracool components is con- 
stantly growing as more searches are performed, they are 
mostly reported in an individual fashion, rather than in 
complete samples that would allow us to asses their true 
frequency. 

Despite the effort that has been devoted to modeling 
over the past decade, evolutionary tracks and atmosphere 
models for ultracool objects still suffer from uncertain- 
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ties. Widely separated L and T dwarf companions to 
nearby stars are also valuable for studying the evolution 
and atmospheres of brown dwarfs. The possibility for 
independent determination of age, mass, and metallicity 
makes these brown dwarfs benchmark objects for cali- 
bration of evolutionary models and atmospheric studies 



lation of 21V1AHS and WISE photometry for 



Kirkpatrick et al. 



(Pinfield et al. 2006 Burgasser et al. 2010a Loutrel et al. 
2011) 



The transition between the spectral types L and T oc- 
curs over a small temperature range of 200-300 K ([Kirk- 
patrick et al.|2000|[Dahn et al.|2002||Tmney et al.|2003| ), 



and is believed to be caused by the depletion ot conden- 
sate clouds, where the driving mechanism for the deple- 
tion is inadequately explained by current cloud models 
(jAckerma n fc Marley|2001||Burgasser et al.|2002||Knapp 



et al.||2004p . I'he complex dynamic behavior ot conden- 
sate clouds of low temperature atmospheres at the L/T 



astrophysics today ( 


Tinney & Tolley 1999[ Bailer- Jones| 


& Mundt 


1999; Gold 


man 2005 Artigau et al.|2U09l|Radi- 


gan et al. 


2012 


)• 



at the L/T transition, and two nearby wide binaries con- 
sisting of M-dwarf prim aries with early L-type and late 
T-type companions FM 

2. SEARCH FOR WIDE BINARIES CONTAINING 
ULTRACOOL DWARFS 

We performed a cross-match between the Two Micron 
All Sky Survey (2MASS) catalog and the Wide-field In- 
frared Survey Explorer (WISE) Preliminary Data Re- 
lease catalog, in order to search for common proper mo- 
tion pairs containing at least one very-low-m ass com- 
ponen t . Following the successful approach of Radigan 



et al. (2008 2009), the correlation of catalogs, calcula- 
tion of proper motions, and identification of co-moving 
stars was done in overlapping sections of 4 deg 2 of the sky 
at a time. For every WISE source, we found the closest 
2MASS match and computed proper motion vectors with 
uncertainties. We then selected stars that had moved at 
> 3<T level compared to all others within the search area. 
Stars within 120" of one another with proper motion 
amplitudes agreeing within 2cr and proper motion com- 
ponents agreeing within la in either right ascension or 
declination were flagged as potential binaries. Towards 
the end of writing this paper, the new all-sky release of 
the WISE catalog became available. We have repeated 
the proper motion measurements for the 4 deg 2 fields 
around the three objects reported here. The main rea- 
son for repeating these measurements is the improved 
astromctry for the faint sources in the new catalog. The 
uncertainties in declination in the WISE preliminary cat- 
alog were computed by adding 0.5" in quadrature to the 
extraction measurement uncertainty to reflect the im- 
pact of the declination bias error known to affect a large 
fraction of sources fainter than Wl ~ 13 mag. This has 
been corrected in the new release. Thus, while the main 
search for binaries has been performed using the prelim- 
inary version, listed proper motions and photometry are 
based on the all-sky release of the WISE catalog. 

11 One of the two binar i es, 2M 1348-1344, was independently 
discovered by |Deacon et al.| | |2012[ l 



recently published a compi- 
350 early 

M- to late T-dwarfs. Mid-infrared WISE photometry 
is particularly suitable for L/T dwarf selection, because 
of their red Wl — W2 colors, which are almost unique 
among astrophysical sources. Additionally, W2 — W3 
color can be used to distinguish between BDs on the one 
side, a nd AGB stars and dus t-obscured galaxies on the 

In the following, we de- 



other (|Eisenhardt et al 
scribe 



|2010p . 

;he procedure for selection of the candidate binary 
systems suitable for spectroscopic follow-up. First, we re- 
quired A, B or C photometric quality flags for J, H, K, 
Wl, W2 and W3 photometry for all the components of 
the potential binary systems. After this initial flagging, 
we were left with ~ 11000 objects. We then rejected all 
candidates for which at least one of the components does 
not satisfy the following conditions: located to the left of 
the solid line in the W2 — W3 vs. W1 — W2 diagram (Fig- 
urefllleft), J - H between -1 and 2, W2 - W3 > -0.5, 
J -W2 > 0.5, and Wl — W2 > 0.35. The latter condi- 
tion was applied to discard a large group of contaminants 
that occupy the color space with Wl — W2 < 0.3 and 
J — W2 < 6 in the right panel of Figure [T] Relaxing this 
condition to Wl — W2 > results in ~ 1200 candidates. 
By applying the cut at Wl — W2 — 0.35, we were left 
with only 19 candidates. Clearly, this restriction makes 
our search insensitive to essentially all M-M pairs, and a 
large portion of L-type companions earlier than L6 will 
be discarded as well. After applying the color cuts on 
either of the components of a candidate binary system, 
we inspected the remaining 19 candidates visually, in or- 
der to remove false positives. In most cases these contain 
blends of unresolved source pairs in WISE, due to signif- 
icantly lower resolution when compared with 2MASS. 

Two color-color diagrams used for the selection are 
shown in Figure [T] Color-code d dots show the colors o f 
M, L and T dwarfs taken from Kirkpatrick et al. (2011 ), 
while the black symbols show the objects reported in this 
paper. In the following subsections we describe the three 
discoveries reported in this paper, along with the known 
parameters obtained from the literature. 

2MASS and WISE photometry for the objects reported 
in this paper are listed in Table [TJ Table [2] contains 
the proper motions and other physical properties of the 
objects discussed in the next sections. 

2.1. 2M06U+3950 

2M06143818+3950357 (hereafter 2M0614+3950) 
was at first identified as a companion to 
2M06143791+3951202. With W2 - W3 « 3, 
Wl - W2 « 0, and J — W2 « 0.75, the latter 
object is marginally consistent with being an early 
M-dwarf. However, an M0 dwarf as faint as J — 16.3 
would have to be at a distance larger than lkpc. By 
checking the optical catalogs of the region, we found 
that the optical source matching 2M06143791+3951202 
has essentially the same position as the 2MASS source 
(offset of - 0.2" from GSC2.2 and USNO-B1 catalog 
positions), whereas at the measured proper motion we 
would expect an offset of ~ 2.5" between the GSC2.2 
and 2MASS, and > 5" between USNO-B1 and 2MASS. 
Furthermore, in the optical images the position reported 
by WISE appears to be located between the source 
matching the 2MASS position and another faint (R~ 19) 
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Fig. 1. — 2MASS and WISE color-color diagrams used for color selection of the MLT dwarf candidates. Diamonds and circles show the 
components of the two new wide binaries 2M0525-7425 and 2M1348 -1344, respectively , while the filled square shows the new L9 dwarf 
2M0614+3950. Dots show the colors of the known MLT dwarfs from [Kirkpatrick et al.| ( |2011| l. The solid line in the left panel is used to 
eliminate extragalactic sources to the right of the line from the bulk of tne M, L,, and T dwarts to the left ( |Kirkpatrick et al.|201l| . We do 
not plot error-bars in cases where they are comparable to or smaller than the plotting symbol. 



source not detected by 2MASS. We thus conclude that 
the WISE source identified here as a potential primary of 
a binary system is most probably a blend. The suspected 
secondary, however, has colors that are consistent with 
a late-L or an early-T ultracool dwarf, and was therefore 
included in our spectroscopic follow-up. 

2.2. 2M0525-7425 

A common proper motion pair here denoted as 
2M0525-7425 consists of the primary 2M05254550- 
7425263 and the second ary 2M05253876-742600 8. Both 
objects were reported by Kirkpa trick et al.] ( |2010[ ) as part 
of a large-area proper motion survey. The secondary 
was observed with CMOS/Gemini during their spectro- 
scopic follow-up, and was assigned the opti cal spectral 
type L2. Th e primary was also reported by |Subasavage| 
et al. (2005) as a high-proper motion source at a photo- 
metric distance of 28.7 pc. Despite the matching proper 
motions, the pair was never reported as a wide binary 
candidate; the primary lacked spectral classification. 

2.3. 2M1348-13U 

A common proper motion pair here denoted as 
2M1348-1344 consists of the primary 2M13480721- 
1344321 (also known as LP 738-14, NLTT 35266, LHS 



2803) and the secondary 2M13480290-1344071. |Reid 
et al. ( 2003 ) published photometry of the p rimary, from 
which they estimate distance of 20.9 pc. |Casagrande| 
et al. ( 2008 ) , again based on photometry, estimated an ef- 
fective temperature of T e g = 2942 ± 58 K. The secondary 
has colors consistent with late T-dwarfs (Figure [T]). 

3. SPECTROSCOPIC FOLLOW-UP 

3.1. EFOSC2/NTT 

We obtained optical spectroscopy of 2M0525-7425A 
(December 23 2011; program ID 184.C-1143(E)) and 
2M1348-1344A (February 29 2012; 088.D-0573(A)) 



with the ESQ Faint Object Sp ectrograph and Camera 
(EFOSC2; |Buzzoni et aT^[l984| on the ESO New Tech- 
nology Telescope (1NTT). Observations were taken in 
long-slit mode with a 1" slit and grism #05, covering 
wavelengths 5200A-9350A with a resolution of 16. 6A. 
The data were flat-fielded, extracted, wavelength- and 
flux-calibrated using the EFOSC2 pipeline provided by 
ESO. The EFOSC2 grism #05 used in this work in- 
troduces stong fringing redwards of about ^7200A. In 
the case of 2M 1348- 1344 A the afternoon dome flats ad- 
equatly removed the fringing, while the effect remains 
present in the spectrum of 2M0525-7425A. The imper- 
fect flat-fielding is due to telescope/instrument flexure, 
which may cause a sligt positional mismatch between the 
fringes in the dome flats and science images. The signal- 
to-noise ratio, as measured by the iraf task splot varies 
between 5 and 18 for 2M0525-7425A, and 5 and 23 for 
2M1348-1344A. 



We 



3.2. SpeX/IRTF 
obtained near-infrared spectrosco py of 



et al. 


2003) 


(lKL'b 


') on 



xna Kea. We used the low-resolution 
prism mode of SpeX to provide a wavelength coverage 
of 0.8-2.5 microns in a single order and an average 
resolution of 120 with a 0'.'5-wide slit. The spectrum 
was acquired using the 0.5- wide slit aligned at parallactic 
angle, with an airmass that ranged from 1.08 to 1.17 
during the observations. The target was observed with 
multiple exposures, dithering in an ABBA pattern along 
the slit, with a total integration time of 2880 sec. After 
the target observation, a nearby AO star at a similar 
airmass to the target was observed for both telluric 
and flux calibration. Internal flat field and argon arc 
lamp calibration frames were obtained at the target 
position for pixel response and wavelength calibration. 
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TABLE 1 
2 MASS and WISE Photometry 



name 




Q 


a 




S 






.1 




H 


K 


Wl 


W2 


W3 


W4 


2M1348-1344A 


13 


18 


07.22 


-13 


11 


32.1 


10.41 


± 


0.02 


9.94 ± 0.02 


9.66 ±0.02 


9.46 ± 0.02 


9.24 ±0.02 


9.11 ± 0.03 


8.72 ±0.33 


2M1348-1344B 


13 


48 


02.91 


-13 


11 


07.2 


16.48 


± 


0.12 


16.09 ±0.17 


> 16.45 


16.15 ±0.07 


14.18 ± 0.05 


> 12.14 


> 9.14 


2M0525-7425A 


05 


25 


45.50 


-74 


25 


26.3 


10.03 


± 


0.02 


9.42 ±0.02 


9.21 ± 0.02 


9.08 ± 0.02 


8.97 ± 0.02 


8.82 ±0.02 


8.53 ±0.15 


2M0525-7425B 


05 


25 


38.76 


-74 


26 


00.8 


15.71 


± 


0.07 


14.97 ±0.10 


14.43 ±0.10 


14.02 ± 0.03 


13.57 ± 0.03 


11.97 ± 0.12 


> 9.46 


2M0614+3950 


06 


14 


38.18 


+39 


50 


35.7 


16.59 


± 


0.14 


15.60 ±0.12 


15.02 ±0.12 


14.33 ±0.03 


13.66 ±0.04 


12.19 ±0.34 


> 9.22 



a coordinates from 2MASS 




^iJAITl) 



Fig. 2.— SpeX spectrum of 2M0614+3950 shown in black, with 
the standard spec tral sequence of ultracool dwarfs at the L/T tran - 
sition (Burgasscr ct al. 2006a Burgasscr 2007] |Looper et al.|2007| ). 
All spectra were normalized at 1.6/um. 

The spectral data was extracted, flat-fielded, wavelength 
calibrated, telluric corrected, and flux-calibrated with 



the Spextool redu ction package ( |Cushing et al.| |2004| 
Vacca et al. 2003). The signal-to-noise ratio ot the final 
spectrum varies between 6 and 36. 

3.3. FIRE/Magellan 

We obtained near-infrared spectroscopy of 2M1348- 
1344B on March 20 2012 with the Folded- port InfraRed 
Echelette (FIRE; |Simcoe et al.|2008"}|2010| spectrograph 
on the 6.5-m Magellan Baade Telescope. In its high 
throughput prism mode, FIRE covers the wavelength 
range from 0.8 to 2.5 /an at a resolution ranging from 
R=500 at J-band to R=300 at Jf-band for a slit width 
of 0.6". A series of ABBA nod pairs taken with expo- 
sure times of 120s per position were used giving a to- 
tal exposure time of 480s. The spectrograph detector 
was read-out using the 4-amplifier mode at high gain 
(1.2 counts per e~) with the SUTR sampling mode. 
We also obtained exposures of a variable voltage quartz 
lamp for flat-fielding purposes and neon/argon arc lamps 
were used for wavelength calibration. The AO dwarf star 
HIP70419 was used for telluric and flux calibration. Data 
were reduced with the new FIRE pipeline tools imple- 
mented in IDL and written by R. Simcoe, J. Bochanski 
and M. Matejek. The signal-to-noise ratio of the final 
spectrum varies between 4 and 14. 

4. RESULTS AND DISCUSSION 

4.1. 2M0614+3950: New brown dwarf at L/T transition 

The spectrum of 2M06 14+3950 is well matched by the 
L9 near-infrared spectral template. In Figure [2] we show 
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Fig. 3. — Color-magnitude di agrams showing ultra cool dwarfs 
with parallaxes from Table 10 of |Dupuy fe Liu| ( |2012| > with errors 
< 0.5 mag. Objects at the L/T transition are shown as purple cir- 
cles. The absolute magnitude of 2M0614+3950 (L9) was calculated 
assuming a distance of d=26.0 ± 1.8 pc; for 2M1348-1344B (T5.5) 
we use d=20.7±1.4pc, and for 2M0525-7425 (L2) d=46.0 ± 3.0 pc. 

the comparison of the spectrum of our ta rget (black) and 
the SpeX NIR standards from L8 to TO ( Burgasser et al. 
|2006a[ iBu rgasser 2007; |Looper et al.|[2007| r~ 

|Dupuy & Liu| ( |2012 ) provide a comprehensive update 
to the absolute magnitudes of ultracool dwarfs with mea- 
sured parallaxes as a function of spectral type, for a wide 
variety of broad-band filters including 2MASS and WISE 
passbands. We use the relations between absolute J, H, 
K, Wl, and W2 magnitudes and spectral type to calcu- 
late the distance to 2M0614+3950. The distance is cal- 
culated as a weighted average of the five estimates, with 
the weights equal to the inverse square of the individual 
uncertainties. These i nclude the rms uncer tainty of the 
polynomial fits given in |Dupuy fc Liu (2012), uncertainty 
of one spectral subtype, and the photometric errors. For 
an object of spectral type L9 we obtain d=26.0± 1.8 pc. 
However, based on the recent analysis of the para llaxes 
for 70 nearby ultracool dwarfs, Faherty et al. ( 2012 ) point 
out that the brightening across the L/T transition is 
more pronounced in reality than when represented by 
the polynomial fits of absolute magnitudes versus spec- 
tral type. Therefore, in addition to the polynomial fits, 
they calculate linear piecewise functions for three spec- 
tral type regi ons. The distance ob tained from the poly- 
nomial fits in Faherty et al. (2012) is in agreement with 
the above estimate. The piecewise fits, however, locate 
the object somewhat closer, at ~ 22 pc (average from 
JHK estimates). We note that distances for objects at 
the L/T transition are uncertain as we are still investi- 
gating their significant scatter in luminosity. 

Figure [3] shows the position of 2M0614+3950 in NIR 



■5 



0.6 
0.4- 
„ 0.2 

3 o.o 

o 

is. 

J -0.2 
-0.4 
-0.6 



2M0525AB 




-0.6 



-0.4 -0.2 0.0 0.2 0.4 
H„cos8 (arcsec/yr) 



0.6 



Fig. 4. — Proper motions of the sources within 1 deg 2 around 
2M0525-7425 and 2M1348-1344. The components of the wide bi- 
nary candidates 2M0525-7425 and 2M1348-1344 are marked with 
the open (primary) and filled (secondary) circles with correspond- 
ing errorbars. 

color-magnitude diagrams, alongside ultrac ool dwarfs 
with para llaxes and errors < 0.15mag from Dupuy fc 
Liu| ( |2012| . 2M0614+3950 is located in the region mainly 
populated by ultracool dwarfs of spectral type between 
~L7 and ~T2. Since the distance to the object was cal- 
culated assuming the absolute magnitude of an L9 dwarf, 
its placement in the correct magnitude range is assured. 
The colors of 2M06 14+3950 are consistent with the L7 - 
T2 spectral range and do not show any peculiarities. 
Using the re lations between effect ive temperature and 

and absolute bolo 



2009) 



spectral type ( Stephens et al.j ^., ...... 

metric magnitude and spectral type ( |Burgasser 2007), we 
obtain T cS w 1350±110K, and log(L/L )=-4 61±0.10. 
The uncertainties are derived by adding in quadrature 
the rms error of the polynomial fits and the error result- 
ing from the uncertainty of one spectral subtype. Com- 



pared to the DUSTY theoretic al isochrones (Chabrier 



et al. 2000 Baraffe et al. 2002), this would correspond 



to an object with a mass of ~ 0.04 — 0.07 Mq, for the 
ages of 1-10 Gyr. A summary of the observational and 
physical properties of 2M0614+3950 is given in Table [2j 

4.2. 2M0525-7425: New wide M/L binary 

Figure [4] shows proper motions of the sources within 
1 deg 2 around 2M0525-7425 and 2M1348-1344. The 
matching radius was set to 2" in order to discard mis- 
matches caused by faint 2MASS sources not detected in 
WISE. The components of the wide binary 2M0525-7425 
are marked with open (primary) and filled (secondary) 
circles with errorbars, which are calculated by combin- 
ing the astrometric uncertainties of the sources in 2MASS 
and WISE catalogs. 

In Figure [5] we show the comparison of the spec- 
trum of the primary (black) and CTIO-4m and CTIO- 
1.5m spectra of dwarfs with spectral types from M2 to 
M4 (purple). The comparison spectra were taken from 
www.dwarfarchives.org, and were chosen to match the 
EFOSC2 wavelength coverage. The spectrum of 2M0525- 
7425A is well matched by the spectra of M2.5 - M3.5 
d warfs. We therefore a ssign it a spectral type M3 ± 0.5. 

|Reid k, Cruz (2002) computed color-magnitude rela- 
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Fig. 5.— EFOSC2 spectrum of 2M0525-7425A shown in 
black, with the spectral sequence of early M-dwarfs from 
dwarf archives . org ovcrplotted in purple. The comparison spectra 
were chosen to have broad wavelength coverage similar to our data, 
and were smoothed to match the resolution of the object spectrum. 
All spectra were normalized at 7500 A. 



tions in the (M v , V-K), (M v , V-I), and (Mj, I — J) 
planes for nearby stars with well-determined trigonomct 
ric parallaxes, which allows us to determine a photo 
metric parallax for 2M0525-7425A. We use JK photome 



try from 2 MASS, V= 
et al. 120051), and 1= 



13.57 (NOMAD catalog; 



Zacharias 



2005}, and 1=11.3 (DENIS cat alog). Combining th e 
results from the three relations in Reid & Cruz (2002), 
we obtain a distance of d=43.9 ± 4.3 pc. The distance is 
calculated as a weighted average of the three estimates, 
with the weights equal to the inverse square of the in- 
dividual uncertainties resulting fr om the scatter of th e 
color-magnitude relations given in Reid & Cruz |2002). 
The photometric uncertainties are either negligible com- 
pared to the fit errors (IJHK), or not available (V), and 
therefore were not included in the calculation. 

The distance to 2M0525-7425B can be estimat ed i n 
the same way as for 2M0614+3950 (see Section |4.1 



The seconda ry was classified as an L2 dwarf ( Kirkpatric, 

distance 



al. 



2010 ): we estimate that it lies at - a 
48.d ± 4.3 pc. The excellent agreement between in- 
dependent distance estimates for the two objects gives 
support to the wide-binary nature of the pair. In the fol- 
lowing, we adopt the weighted average d=46.0±3.0pc as 
the distance to the 2M0525-7425 binary, with the weights 
equal to the inverse square of the component distance un- 
certainties. At this distance, the projected separation of 
43.9" ± 0.1" is equivalent to 2020 ± 130 AU. The uncer- 
tainty of the projected separation combines the astro- 
metric errors of the binary components from 2MASS. 

Using the re lations between effect ive temperature and 
spectral type (Stephens et al. 2009), and absolute bolo- 



G 



Muzic et al. 



metric magnitude and spectral type (Burgasser 2007), we 
obtain T eS w 1970 ± 120 K, and log(L/L )=-3.79±O.ll 
for the L2 comp anion. The spectral type f or this source 

, which gives 
therefore as- 



is adopted from Kirkpatrick et al 
no formal errors on spectral types 



(2010) 
We 



sume an uncertainty of 0.5 subtypes, that is equal to 
the step in their spectral typing scheme. The resulting 
uncertainty was added in quadrature to the scatter of 
the polynomial fits. For the primary, we use relations 
between effective tempe rature and color determined by 
Casagrande et al. (20081) for a large sample of nearby M 
dwarfs with accurate optical and infrared photometry. 
We calculate T c g as the average of the estimates based 
on six optical-to-near- infrared colors (V and R versus J, 
H, and K), with the uncertainty equal to the standard 
deviation of these six estimates. The R magnitude orig- 
inates from the HST Guide Star Catalog-II (GSC-II)Fj 
We obtain T eff « 3320 ± 50 K. From the relation be- 
tween the 2MASS and bolometric magnitudes we obtain 
log(L/L )=— 1.60 ± 0.03. The uncertainty takes into ac- 
count the scatter of the linear relations between JHK 



and bolometric magnitudes (Casagrande et al. 2008), 



p hotometric errors, an d the distance uncertainty. 



Lepine et al 
index Q using t 



(20071 defined the metallicity dependent 
re CaH2, CaH3 and Ti05 molecular band 



heads in th e optical, recently re-calibrated by [Dhital 



et al.j (|2012J). From the EFOSC2 spectrum of 2M0525- 
7425 A we find Cl=0.94 and Cd=0.89, placing the object 
among the dwarf metallicity class. However, while the 
(-index can be very useful for discriminating between 
different metallicity classes, its correlation with [Fe/H] 
within the dwarf metallicity class (C > 0.825) is weak and 
shows a significant scatter ( Lepine et al.|2012 ). The loca- 
tion of M dwarfs in the (V — K$) - Mk s color-magnitude 
diagram has been s hown to correlate with metallicity 
dBonfils et al. || 2005 | jJohnson fc Apps||2009| I Schlaufmanl 
&; Lau ghlin 20 10[). Using the relation from pchlautman] 
Hz haughiinf (|2010[ ), who improved upon previous cali- 
brations, and adopting the average distance to the bi- 
nary d=46.0±3.0pc, we obtain [Fe/H]=-0.01 ±0.03 for 
2M0525-7425A. The final uncertainty combines the un- 
certainties of the distance and the Kg magnitude, while 
the error of the V magnitude is not available. Note that 
Neves et al.|(|2012[) find the rms dispersion of the Schlauf- 



man fc La ughlin (2010) relation to be 0.19 ± 0.03 dex 

| West et al.| ( |2008[ ) showed that the decrease in the 
fraction of magnetically active stars (as traced by Ha 
emission) as a function of the vertical distance from the 
Galactic plane can be explained by thin-disk dynamical 
heating and a rapid decrease in magnetic activity. The 
timescale for this rapid activity decrease changes accord- 
ing to the spectral type. From the lack of Ha emission 
in the spectrum of 2M0525-7425A, we place a lower limit 
for its age at 1.5 Gyr. 

Assuming an age of 1-10 Gyr, the component masses 
can be estimated from evolutionary models. By com- 
paring the derived absolute JH Ks magnitudes wi th the 
BCAH98 evolutionary models (Baraffe et al. 19981, we 
find a mass of 0.46 ± 0.01 M for the primary. The 
adopted uncertainty was determined by varying the ab- 
solute magnitudes within their corresponding uncertain- 

12 Photographic R f magnitude from GSC-II equival ent to the 
Johnson-Cousins R c | |Reid et al.|1991| |Kirk et al.|2009) . 



ties. This result is in agreement with the value 0.45 ± 
0.01 M obtai ned using the empiric al mass-luminosity re- 
lationship of Delfosse et al. (2000). For the secondary, 



we find a mass of 0.070^q'qiq M , with an error reflect- 
ing the model grid spacing. The mass ratio of 2M0525- 
7425AB is 0.15 ± 0.02, which is relatively low for a low- 
mass pair. The mass ratios are discussed in Section [B~2 
A summary of observational and physical properties o 
2M0525-7425AB is given in Table |2j 

4.3. 2M1348-1344: New wide M/T binary 
In the left panel of Figure [6] we show the compari- 



son of the EFOSC2 spectrum ofthe primary (black) and 
the CTIO-4m and CTIO-1.5m spectra of dwarfs with 
spectral types from M4 to M5 (purple). The compari- 
son spectra were taken from www.dwarfarchives.org, and 
were chosen to match the EFOSC2 wavelength coverage. 
Based on this plot, we assign 2M1348-1344A a spectral 
type of M4. 5 ±0.5. The FIRE spectrum ofthe secondary 
is shown in black in the right panel of Figure |6l along 
with the SpeX NIR standards from T4 to T8 (Burgasser 



et al.||2004 
provide goo 



2006a]. Spectral templates T5 and TB both 
iTTEo the spectrum of 2M1348-1 344B. As a 
check , we ca lculate the spectral indices from |Burgasser| 

of 



^tEe 



major 



et al.| ( |2006a|) , that are based on the strength < 
H2O and CH4 bands in the near- infrared. Five indices 
yield spectral types from T5 to T7, with the average of 
T5.7. We type 2M1348-1344B as T5.5 ± 1.0. 

The distance to the primary can be ca lcula ted in the 
same way as for 2M0 525-7425A (Sectio n |4.2[ ) , adopting 
V=15.1 and 1=12.1 ( |Reid et al.|[2003| ). This gives us 
d = 21.5 ± 2.1 pc. The distance to the T-dwarf was cal- 
culated in the same way as for the L2 and L9 dwarfs 
from the previous two sections, based on J, H, Wl and 
W2 magnitudes and by adopting the spectral type T5.5. 
We obtain d= 20.0 ± 2.0 pc . The excellent agreement 
between the independent distance estimates for the two 
objects supports a binary classification. In the follow- 
ing we adopt the weighted average d=20.7 ± 1.4 pc as 
the distance to 2M1348-1344 binary. At this distance, 
the projected separation of 67.6" ± 0.3" is equivalent to 
1400±90AU. 

Following the same procedure as in Section [4~2| we ob- 
tain T cS w 2960 ± 20 K, and log(L/L )=-2.42 ± 0.03 
for the primary. The activity lifetime for M4 dwarfs is 
4.5 ± 0.5 Gyr. From the lack of Ha emission in the spec- 
trum of 2M1348-1344A, we place a lower limit for its 
age at 4 Gyr. The (-index of 2M 1348- 1344 A is found 
to be very close to unity ((lrs = 1-08, Cd = 1-06), 
placing the object into the dwarf metallicity clas s . The 
photometric relation of Schlaufman & Laughlin ( 2010[ ) 
results in [Fe/H]=— 0.33 ± 0.08. Here we have to take 
into account the rms dispersion o f the relation tha t 
equals 0.19 ± 0.03 , as calculated by Neves et al. (2012). 
Our result for the metallicity of 2M1348-1344A is there- 
fore in agreement with the values found for other M- 
dwarfs in the solar n eighbourhood, which haye the mean 
[Fe/H]=-0.17±0.07 ( [Schlaufman fc Laughlin|2010| ). The 
estimate for the mass of the primary gives 0.18±0.01 M 
(BCAH98, age 5-10 Gyr). The s ame number is obtaine d 
by following the relationship of Delfosse et al. ( |2000[) . 
For the T5.5 secondary, we obtain T eS sa 1070 ± 130 K 
and log(L/L )=-5.O3 ± 0.21. Using COND03 ((Baraffe 
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Fig. 6. — Left: EFOSC2 spectrum of 2M1348-1344A shown in black, with the^M^dwarfspectral sequence from dwarfarchiveSjOrg 
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match the resolution of the object spectrum. All spectra were normalized at 7500 A. Right: FIRE spectrum of 2M1348-1344B shown 
black, with the standard spectral sequence of T dwarfs (Burgasscr ct al. 2004 2006a). The higher-resolution F1KE spectrum was smoothed 
fc^naTcn^reTescTniUoirTFtEe^specCT 



et al.|2003 ) evolutionary models, we find the mass of the 
secondary to be 0.04 ± 0.01 M for ages 5-10 Gyr. As in 
the case of 2M0525-7425, we find a relatively low mass 
ratio of 0.22 ±0.06. A summary of the observational and 
physical properties of 2M1348-1344 is given in Table [2] 

2M1348-1344B is one of only a handful of known late 
T (>T5) wide companions to low mass stars, and one of 
only six late T dwarfs separated by more than 1000 AU 
from their primary stars. The other objects include Gl 
570D (K4 /M1.5/M3/T7 system K4 and T7 separated by 



3600 AU; [Burgasser et al.pOOO ), e Indi Bb (K5/T1/T6 
system, K5 separated b y~1460 AU from the T1-T6 clos e 
binary; iScholz et all 120031 iMcCaughrean et al.||2004h, 
Ross 458TTPV10/M7/T8.5, MO an d T8. 5 separated 1170 
AU; |Scholz||20T0l |Goldman et al.| |20l0 Burg asser etH] 



2010b Burningham et al. 



IS 



with separation ot 1200 A 
2011 ), and G204-39 B (M3/T6.5 separate^ 
Faherty et aT|2010| ). 



H ip737 86B (K5/T6p 
olz||2010 ™ 




Murray et a 
by 2685 AU; 



5. DISCUSSION 

5.1. Stability of the wide binaries 

Figure [71 shows the total mass versus separation for a 
sample ofbinary systems. We show the suggested em- 
pirical li mits for the stabilit y of VLM multiples from th e 



studies of |Reid et al. (20011 and Burgasser et al. (20031, 
based on objects that were known at the time of the re- 
spective works. Neither cutoff seems appropriate for the 
collection of widely separated systems plotted (see dis- 



10 3 10* 10 1 

Separation UCD to Star (AU) 

Fig. 7. — A plot of the total mass vs. binary separation. Ob- 
jects marked with filled circles are tight low-mass systems. Wide 
systems (separation >100 AU) containing a UCD companion are 
marked as five point stars and those marked as blue squares are 
systems containing a tight or widely separated UCD with an age 
< 500 Myr. Objects marked by open circles come from stellar 
companion catalogs. The new wide systems reported in this work 
are marked as red five-point stars. We have over-plotted the |Reid| 
|et al.| ([2001 ) curve (center) which distingui shed the cut-off for the 
tormation of wide stellar pairs as well as the Burgassc r et al.| | |2003[l 
line (specific for Mtot <0.2 Mq field systems) and t he log-normal 
limitation for systems with Mtot > 0.3 Mw found in[Dhital et al.| 
| |2010[ l. Data were taken from Figure 9 in Fa herty et al.| l |2Ull[ |, see 
rctcrcnccs therein for more details. 



Faherty et al. 


2010 


and 


Dhital et al. 


2010) 



cussions m 

The log-normal limitati on on the separation of binaries 



et al. (2010) used the Jeans length instead of the fiducial 
value of 300 AU to repeat the calculation for two extreme 
mass ratio cases. Both 2M0525-7425AB and 2M1348- 



found by |Dhital et al. (2010) does encompa ss 2M0525- 
7425AB and a ll low-mass objects in Figure[7J Zuckerman| 
fc Song| ( |2009[ ) derived a cutoff for the binding energy of 
a system tormed by fragmentation as a function of the to- 
tal system mass, assuming separations of 300 AU. Since 
a number of systems found in the field and yo ung clus- 
ters are found at separations larger than 300 AU, |Faherty| 



Faherty et al. 



Dhital et a 



1344 are found abov e the Zuckerman & Song ( 2009 1 and 



(2010) limits. 
12010) calculate the average binary life- 



times, that generally depend on the total mass and sep- 
aration (see their Figure 19). With the parameters cal- 
culated in the previous sections, both systems reported 
here appear to be stable enough to survive longer than 



Muzic et al. 



TABLE 2 

Physical properties of 2M1348-1344AB, 2M0525-7425AB, and 2M0614+3950 



2M1348-1344A 



2M1348-1344B 



2M0525-7425A 



2M0525-7425B 



2M0614+3950 



2MASS designation 
fi a cosS (mas yr ) 
US (mas yr -1 ) 
Spectral type 
Distance (pc) 

Teff(K) 

Mass (M Q ) 
Age (Gyr) 
[Fe/H] d 

Angular separation (") 
Projected separation (AU) 
Binding energy (10 41 erg) 
Mass ratio 



13480721-1344321 13480290-1344071 05254550-7425263 05253876-7426008 06143818+3950357 



-680 ± 28 
-511 ± 32 
M4.5 ± 0.5 b 
21.5 + 2.1 
2960 ± 20 
0.18 + 0.02 

> 4.0 
-0.33 + 0.08 

67.6 + 0.3 
1400 ± 90 
0.91 + 0.26 
0.22 + 0.06 



-657+ 23 
-535 + 33 
T5.5±l c 
20.0 + 2.0 
1070 ± 130 
0.04 + 0.01 



270 ± 18 
308 ± 23 
M3.0 ± 0.5 b 
43.9 + 4.3 
3320 ± 50 
0.46 ± 0.01 

> 1.5 
-0.01 ± 0.03 

43.9 + 0.1 
2020 ± 130 
2.81 + 0.45 
0.15 + 0.02 



279 ± 18 
310 + 23 

L2 ab 
48.0 + 4.3 
1970 ± 120 



070+ 005 
u.u(u_ 010 



-46 + 21 
-265 ± 23 
L9 + l c 

26.0+ 1.8 
1350+ 110 
0.04 - 0.07 



a |Kirkpatrick et al.| l |2010( 
b optical spectral type 
c NIR spectral type 

Photometric relation from |Schlaufman & Laughlin| l |201o| 

10 Gyrs. Based on the component proper motions, ex- 
cellent agreement in the independent distance estimates 
of each component, and stability arguments, we conclude 
that 2M0525-7425AB and 2M 1348-1 344 AB are genuine 
binary systems. 

5.2. Mass-ratios 

Both wide binary systems reported in this work have 
relatively low mass ratios(g = 0.22 for 2M1348-1344 and 
q = 0.15 for 2M0525-7425). By analyzing a catalog of 
1342 very- wide (projected separations of > 500 AU) low- 
mass ( at least one m i d-K to mid-M dwarf component) bi- 



naries, 
ratios 



Dhit al et al.| (2010) find the distribution of mass 
o be strongly skewed toward equal-mass pairs: 
85.5% of pairs have masses within 50% of each other. 
The authors conclude that this does not result from ob- 
servational biases, i.e. the survey is sensitive to pairs 
with much lower mass ratios. A similar result is obtained 



by ( Reid & Gizis 1997 ) in the volume- complete 8-pc sam- 
ple, consnT^TEecIinostly (80%) of M-dwarfs. However, the 
sensitivity to mass ratios in both works is reported for 
the entire samples, rather than for M dwarfs alone, and 
the sensitivity to low mass ratio systems will fall off with 
lower-mas s primaries. In a study of multiplicity among 
M-dwarfs, [Fischer fc Marcy| ( |1992"| ) found that the distri- 
bution of mass ratios is flat and perhaps declines at low 
mass ratios. A number of low-g binary systems contain- 
ing ultracool components are known to date, some with 



the mass-ratios as low as 0.05 (e.g. Faherty et al. 2010 



and references therein). A complete sample sensitive to 
ultracool dwarfs of the very lowest masses is needed to 
assess the distribution of mass ratios in this regime. 

5.3. Higher-order multiplicity 

It has been found that binaries wider than ^100 AU 
have a high rate of tertiaries, both in the stellar and 
substellar regimes, and that the binary fraction of ul- 
tracool companions found in multiple systems ap pears 
to be significantly larger than that of field dwarfs (Bur 



asser et al.|2005||Tokovinin et al.|2006] |Law et a"L]? 



'aherty et al. 2010p . This behavior has also been pre- 
dicted by some dynamical models ( Sterzik fc Durisen 
2003] IDelgado-Donate et al.| [20041 |Whitworth fe Sta- 



matcllos 2 006[ [Kouwenhoven et al. |2010 ). Furthermore 
dwarfs with spectral types of late L or early T show e 



higher mult i plicity fraction compared to other types ( Liu 
et al. 20061 IBurgasser et al . 2006b; but see also |Gola 



man et al.||2008p . This has been attributed to the rapid 
evolution through the L/T transition phase, implying 
fewe r sources per spe ctral subtype for a given field sam- 
ple ( Burgasser||2007 ). Unresolved binaries can be iden- 
tified by an apparent over-luminosity for a given spec- 
tral subtype, whe re an independent d istance estimate 



is available (e.g. Dupuy & Liu 2012) or by spectral 



deviat io ns from template spectr a (e.g. Burgasser et al. 



2010aJ. |Burgasser et al. (|2010aJ identified several com 
mon trends that stand out when comparing the combined 
near- infrared spectra of resolved L/T binaries with the 
template spectra of single objects. 

The SpeX spectrum of the L/T transition dwarf 
2M0614+3950 (Figure [2} is closely matched by the L9 
spectral template, and goes not show any of the devia- 



tions pointed out in |Burgasser et al~ (2010a). For the two 
wide binaries, 2M0525-7425AB and 2M1348-1344AB, the 
distance estimates of the primaries are in excellent agree- 
ment with the independent distance estimates of the 
secondaries. For example, by placing the T5.5 dwarf 
2M1348-1344B at the distance of the primary augmented 
by its ler uncertainty +0.25 in its absolute J-magnitude, 
and -0.15 in its J — H color, still consistent with the 
derived spectral type. Also, all of the component spec- 
tra presented here do not appear to deviate significantly 
from the spectral templates. We note that the classifica- 
tion o f 2M0525-7425B originates from Kirkpa trick et alT] 
(20101, who do not show a spectrum in the paper, but 
neither report any peculiarities. We conclude that, while 
we have no reason to suspect any unresolved binarity 
in the objects discussed here, high resolution imaging 
and/or radial velocity study is necessary to prove this 
claim. 



5.4. Predictions from theory 

Among the multiple systems produced in hy drodynam - 
ical simulations of star cluster formation by 
about 10% are stars (> 0.1 M ) with VLM 



Bate (2009}, 

companions. About 30% of these star/VLM binary sys- 
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terns have semi-major axes larger than 1000 AU, and 
several objects with primary masses similar to those of 
the two systems reported here are created. Although 
there appears to be no statistically significant depen- 
dence of the frequency of star/VLM binaries on primary 
mass, the typical separation of such binaries seems to 
increase strongly with the primary mass. As for the 
mass ratios for the systems with 0.1 - 0.5 M primaries, 
Bate (2009) predicts a slight preference towards near- 
equal mass systems. However, a significant number of 
systems with low-mass ratios is produced in the simu- 
lations, with about 15% of the binaries having q < 0.2. 
Another prediction of this model is a corrrelation be- 
tween mass ratio and separation for the binaries, with 
closer systems having a preference for equal masses. The 
median mass ratio for binary separations in range 1000 - 
10 4 AU is 0.17, comparable to the mass rat ios observed 
in our two binary sy stems. Simulations by |Stamatellos| 
fc Whitworth] ( |2009j ) show that very wide (up to 10 000 
AU), low-mass (secondary mass < 80 Mj) companions 
can be formed in fragmenting stellar disks. In the sim- 
ulation of a 0.7 M star with a disk of an equal mass 
and a radius of 400 AU, wide (400-10 000 AU) binary 
companions to the central star are predominantly brown 
dwarfs. An alternative mechanism to form wide binaries 



was proposed by Kouwenhoven et al. (2010) and Moeckel 
& Clarke (2011). Their IN-body simulations of star clus 



ters in dissolution show that binaries with separations 
larger than 1000 AU can be formed by dynamical cap- 
ture of initiall y unbound objects. While t he wide binary 
fraction from Kouwenhoven et al. (2010) varies signifi- 
cantly with the initial conditions, by integrating over the 
initial cluster mass distribution they predict a wide bi- 
nary fraction of a few percent in the Galactic field. For 
wide systems with se parations > 1000 AU and primary 



masses below 1.5 M , Kouwenhoven et al. (2010) models 



predict the increase in the number ot the systems with 
decreasing the mass ratio. The mass-ratio distributions 
peaks around q = 0.1 — 0.3. 

The existence of binaries with total system masses, 
mass ratios and separations similar to the two systems 
reported here is clearly supported by predictions of dif- 
ferent theoretical models. A proper statistical analysis 
of complete wide binary samples is needed to asses their 
overall frequency and distribution of properties such as 
mass ratio, and to be able to compare these observational 
results with the predictions of the models. 

6. SUMMARY AND CONCLUSIONS 

We have presented the results of our search for com- 
mon proper motion pairs containing ultracool compo- 
nents, based on 2MASS and Preliminary WISE data cat- 
alogs, followed by optical or NIR spectroscopy of individ- 
ual candidate components. We report the discovery of 
two wide binaries containing ultracool components, and 



one ultracool dwarf at the L/T transition. 

2M1348-1344AB consists of a M4.5 primary with a 
T5.5 companion, at a distance ~ 20 pc from the Sun, 
and at separation of 1400 AU. An M3 primary and its L2 
companion comprise the system 2M0525-7425AB, found 
lying at a distance ~ 45 pc, and separated by 2000 AU. 
Based on the matching proper motions and excellent 
agreement in distance estimates of the components, as 
well as on the fulfilled stability criteria, we conclude that 
both systems are genuine low-mass binaries. Both pri- 
maries are found to have metallicities close to that of 
our Sun, and lack Ha emission, allowing us to place the 
lower limits on the ages of the two systems at 1.5 Gyr 
(2M0525-7425) and 4 Gyr (2M 1348- 1344). We also re- 
port the discovery of 2M06 14+3950, an ultracool dwarf 
with spectral type L9, lying at a distance ~ 26 pc. 

The three objects reported here add to the com- 
pendium of ultracool objects discovered by WISE, which 
is certain to grow with the recent survey's all sky release. 
Further characterization of the L/T transition dwarf and 
the components of the two new binary systems will pro- 
vide stringent tests for understanding of ultracool atmo- 
spheres, as well as for the formation models of BDs and 
BD binaries. 2M1348-1344AB is particularly important 
in this sense, being one of a handful of known late-T very 
wide companions to low mass stars. 
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